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12 We do not demand from the author to transfer his own copy- 
right to us. Instead, we only consider papers resulting 
only from original research work, and only if the author 
does assure us that he owns the copyright of his paper and 
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Visual EM Simulator for 3D Antennas: 
VEMSA3D - FLOSS for MS Windows 



C . A . Koutsos, N.I. Yannopoulou, P.E. Zimourtopoulos * 

Antennas Research Group, Slovak Republic [l] 
Antennas Research Group, Austria - Hellas [2, 3] 

Abstract 



This paper introduces the FLOSS Free Libre Open Source 
Software [VEMSA3D] , a contraction of "Visual Electromagnetic 
Simulator for 3D Antennas", which are geometrically modeled, 
either exactly or approximately, as thin wire polygonal stru- 
ctures; presents its GUI Graphical User Interface capabili- 
ties, in interactive mode and/or in handling suitable formed 
antenna data files; demonstrates the effectiveness of its use 
in a number of practical antenna applications, with direct 
comparison to experimental measurements and other freeware 
results; and provides the inexperienced user with a specific 
list of instructions to successfully build the given source 
code by using only freely available IDE Integrated Develop- 
ment Environment tools-including a cross-platform one. The 
unrestricted access to source code, beyond the ability for 
immediate software improvement, offers to independent users 
and volunteer groups an expandable, in any way, visual an- 
tenna simulator, for a genuine research and development work 
in the field of antennas, adaptable to their needs. 

Keywords 



FLOSS, antennas, modeling, simulation 



Introduction 

A lot of amazing visual EM 
software simulators, both com- 
mercial and freeware, exist 
for many years now. However, 
to the best of authors' know- 
ledge, none of these simula- 
tors has been ever released 
under an Open Source license- 
paid, free gratis or free li- 
bre. This situation, as it 
was discussed recently, pre- 



cludes the independent users 
and volunteer groups with 
limited resources, like the 
authors' nonprofit group, from 
the genuine, state-of-the- 
art, scientific research, sin- 
ce, neither the code improve- 
ment, nor its adaptation to 
specific needs, is possible 
[ 1 ], [ 2 ]. 

For that reason, the au- 
thors decided to develop and 
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release under the GNU public 
license their own Visual Elec- 
tromagnetic Simulator for 3D 
Antennas [VEMSA3D] , although 
it certainly is a less well- 
equipped software application 
for the moment-but still a 
fully expandable one. That re- 
lease is now perfectly per- 
missible, because in 2005 
NASA released in the public 
domain the FORTRAN source co- 
de of the well-known MoM Me- 
thod of Moments Thin-Wire Com- 
puter Program, by J. H. Ri- 
chmond [3], [4]. This is ex- 
actly the code on which the 
authors' group based its free- 
ware simulators: [RichWire], 
a CLI Command Line Interface, 
and [DA], a MS Quick-Win FOR- 
TRAN derivative of it. These 
two simulators are under un- 
interrupted revision, impro- 
vement, expansion, and rede- 
velopment, since a long time 
ago [1], [2]. Therefore, the 
[Rich-Wire] FORTRAN code was 
translated-entirely, line-by- 
line, without using any paid 
or free translator-to C++, to 
form the core of scientific 
EM computations in [VEMSA3D] . 

On the other hand, the au- 
thors' group requirements for 
scientific accuracy in visual 
representation of the produ- 
ced EM simulation antenna re- 
sults from its simulators, we- 
re already enforced the soft- 
ware development of the Vir- 
tual Antennas, that is the 
Virtual Antennas laboratories. 



in VRML [5], [6], the FLOSS 

application [RadPat4W] for an- 
tenna radiation pattern pre- 
sentation, in MS VB6 [2], 
[7], as well as, the recently 
developed visual antenna ap- 
plication for the Wolfram De- 
monstrations Project in Ma- 
thematica [8] . The visualiza- 
tion ideas implemented in the 
aforementioned graphics ap- 
plications were also expres- 
sed-from the scratch-in C++, 
to use the Open Source cross- 
platform [wxWidgets] library 
with OpenGL and form the core 
of scientific EM graphics in 
[VEMSA3D] [9] . 

The authors, having taken 
into account that no famil- 
iarization with software use 
is possible without getting 
hands-on experience, restrict 
themselves to a brief discus- 
sion of the current [VEMSA3D] 
characteristics. The code, the 
antenna applications data, as 
well as, any other informa- 
tion, referenced open bugs or 
future code releases, will be 
always available in authors' 
group repositories, at 
"www. antennas . gr /floss" 
and in GoogleCode website, at 
"http : //code . google . comp/p/rg 
a/" . 

GUI Interactive Mode of 

Operation 

It is assumed that the 
user has already some experi- 
ence in the sketching of a 
polygonal wire outline model, 
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for an antenna under consid- 
eration, consisting of num- 
bered wire segments and no- 
des, with their 3D space co- 
ordinates, including the po- 
sitions of any antenna cir- 
cuit elements, that is volt- 
age generators and lumped 
loads, and s/he wants then to 
key-in these model and cir- 
cuit data into [VEMSA3D] us- 
ing its GUI Graphical User 
Interface in interactive mo- 
de . 

The GUI main window is 
shown in Fig. 1, while Fig. 2 
shows the menu items along 
with their available submenu 
options numbered from 1 to 7. 

The GUI itself is divided 
in three panels named [An- 
tenna in space (3D)], [Anenna 
Elements], and [Data-sets]. 



The function of each panel is 
briefly described in the fol- 
lowing . 

The [Antenna in space 
(3D)] viewing panel is used 
to project all the generated 
3D and 2D graphics. By de- 
fault, the application starts 
with a simple linear dipole 
antenna in space, as it is 
shown in Fig. 1 from a view- 
point on the diagonal of tri- 
hedral angle OXYZ, where the 
usual letters associated with 
the Cartesian axes do not ex- 
ist. Instead, a one-to-one 
correspondence implicitly ex- 
ists between (X, Y, Z) axes 
and (R, G, B) colors [5]. The 
3D image of any antenna can 
be manipulated through the 
[Study] options, as shown in 
Fig. 2.3. 




Fig. 1: GUI: The main window 
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Options MIS 
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Simulation Options Help 
• Antenna in space (3D) 4 
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Fig. 2: Main window: The unfolded menu items 



The [Antenna Elements] pa- 
nel contains four tabs, cor- 
responding to model and cir- 
cuit data, in four self-ex- 
plained building blocks. In- 
stead of attempting to de- 
scribe the block data in gen- 
eral the presentation conti- 
nues with the default data: 

[Nodes] : 2, with coordi- 
nates 1: (0.2, 0.1, 0.25), 2: 
(0.2, 0.1, -0.25), in meters 
[m] or in terms of wavelength 
A [wl]. 

[Wires] : 1, starting at 
node: 1, ending at node: 2, 
that is a directed segment. 

[Segments] is either: (a) 
a positive integer, (b) a 
zero or (c) a number between 
0.00 and 0.25, to respec- 
tively define the division of 
this wire in segments: (a) in 
the indicated number of seg- 
ments, (b) in a calculated 
number of segments of length 
no longer than 0.05A (that is 
10 segments, in this case), 



(c) in segments with length 
no more to that number. No- 
tably, no segment can be 
longer than 0.25A [3]. 

[Generators] : 1, connected 
in the wire:l, at a distance 
of 0.5 times this wire 
length, from its starting no- 
de, with rms value of 1^0° 
[V] . There can be only one 
generator on each wire. 

[Loads]: Has the same stru- 
cture as [Generators] . Resis- 
tors, inductors and capaci- 
tors are inputted in [Q], [H] 
and [F], respectively. It is 
empty in this case, since no 
lamped loads exist. 

Notably, the GUI has been 
developed in such a way that 
efficiently supports the in- 
teractive handling of the an- 
tenna modeling by the user on 
the screen. Thus, the selec- 
tion of an element in any of 
these four tabs simultane- 
ously highlights this element 
on the viewing [Antenna in 
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space (3D)] panel-if the [Se- 
lect] option under [Study] 
(Fig. 2.3) has been already 
chosen by the user, and this 
is the default selection in 
any case-and conversely: the 
selection of any element on 
this viewing [Antenna in spa- 
ce (3D)] panel highlights the 
element in its table. Conse- 
quently, a full control on 
the antenna modeling struc- 
ture is achieved. 

The [Datasets] panel may 
contain independent sets of 
model and circuit data for va- 
rious antennas. Such a Data- 
set has its own [Dataset Si- 
mulation Options], which are 
chosen through submenu of 
Fig. 2.6, as the default se- 
lection of them is shown on 
the left part of Fig. 3. 

A [Wire conductivity] of -1 



corresponds to perfect wire 
conductivity, while the [Spa- 
ce conductivity] of 0 and the 
[Space dielectric constant] 
of 1 defines the Free Space 
EM environment. The [Integra- 
tion steps] accepts a posi- 
tive number for the steps of 
approximated integration used 
in the MoM impedance computa- 
tions-a zero means closed- 
form integration. 

Notably, the [Maximum seg- 
ment length] affects the di- 
vision of wires: the smaller 
this number is, the more seg- 
ments will be considered. It 
is currently known that there 
is an implied, still program- 
matically unimproved, priori- 
ty of the number of segments 
defined in [Wires] tab, over 
this selection-that is a bug, 
which may crash the applica- 
tion . 



Dataset Simulation Opti 
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[MHz] 


Wire radius : 


0.001 


[m] 


Wire conductivity : 


1 _1 
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Space conductivity : 


0 


[S/m] 


Space dielectric constant : 


1 


]B 


Maximum segment length : 
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Integration steps : 


0 


]0 


3D angle step : 


5 


[deg.] 


2D angle step : 


1 


[deg.] 



OK ~j | Cancel | | Default ] 
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Simulation started. 
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Running Dataset : Dataset 


20:20:35.687 


SEG started. 
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SEG finished. 


20:20:35.703 


SORT started. 


20:20:35.750 


SORT finished. 


20:20:35.750 


SGANT started. 


20:25:32.718 


SGANT finished. 


20:25:32.718 


GANT1 started. 


21:03:11.906 


GANT1 finished. 


21:03:11.906 


PATTERNS started. 


21:03:22.500 


PATTERNS finished. 


21:03:22.500 


Simulation finished. 


21:03:22.500 


Total time: 00:42:47.047 



Close 



Fig. 3: Windows: [Dataset Simulation Options] - [Status] 
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The EM simulation of an an- 
tenna starts by selecting [Si- 
mulate Current Dataset], un- 
der [Simulation] (Fig. 2.5), 
while its progress is shown 
in the [Status] window, such 
as the one at the right of 
Fig. 3, which concerns the 
default dipole. Notably, mul- 
tiple Datasets can be simu- 
lated, one after the other, 
through the [Simulate Marked 
Datasets] shown in Fig. 2.5. 

The course of the simula- 
tion may be roughly described 
in five steps: 

(1) Wires are automatically 
segmented, using the provided 
parameters, to produce the fi- 
nal model structure of points 
and segments 

(2) Model structure is ana- 
lyzed and adjacent segments 
are combined to form dipole 
current modes with sinusoidal 



distribution 

(3) Impedances are calculated 
and the MoM algebraic system 
of equations is formed 

(4) The system of equations 
is solved and the segment cur- 
rents are calculated 

(5) 3D and 2D radiation pat- 
terns are calculated, as well 
as other useful results like 
those shown in Fig. 4. 

In Fig. 5, the [Viewing 
Options] under the [Options] 
(Fig. 2.6) is given with all 
of the available options for 
the graphics. Through [Graph] 
menu item (Fig. 2.4), a vari- 
ety of 3D and 2D plots, which 
include either the normalized 
electric far field E radia- 
tion pattern or the relative 
radiation intensity U in dB, 
as well as their 0 and cp 
parts, can be illustrated. 




Fig. 4: Window: [Results] - Table: [Power] 
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The [Antenna segments in 
space (3D)] illustrates the 
final antenna structure with 
points, segments and segment 
currents in amplitude and 
phase, as shown in Fig. 6, 
where the resulting current 
distribution on default di- 
pole it seems to be sinu- 
soidal, that is as it was ex- 
pected to be. The well-known 



2D intersections of the 3D 
radiation patterns by the 
three main planes yOz, xOy, 
zOx, as well as these 3D pat- 
terns themselves are also 
shown for the default dipole. 

Notably, modeled antennas 
can be imported and exported 
in [RichWire] data file for- 
mat (Fig 2.1) . 




Fig. 5: [Viewing Options] 




Fig. 6: GUI: Default dipole results 



Build under MS Windows XP 

To build [VESMA3D] MS Win- 
dows executables, the use of 



IDE CodeBlocks with GCC/Mingw 
is suggested [10] : 
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Set below as X:\ the system disk, e.g. X:\ -> C:\ 

1 Download (~20MB) , Install (~110MB) : 

http : //sourceforge . net/pro jects/codeblocks/f iles/Binaries/ 
[View all files] [Binaries] [8.02] 
codeblocks-8 . 02mingw-setup . exe 
Select: [Full], [Run], Compiler: [GCC/Mingw] 

2 Download wxWidgets (~20MB) , Extract (~120MB) : 

http : //sourceforge . net/pro jects/wxwindows/f iles/2 . 8 . 11/ 
[View all files] wxWidgets-2 . 8 . 11 . zip 
wxWidgets-2 . 8 . 11 . zip-> X:\ [Yes to All] 

3 IDE CodeBlocks, Open: 

X : \wxWidgets-2 . 8 . ll\include\wx\msw\setupO . h 
Replace, at line 1006: 

#def ine wxUSE_GLCANVAS 0 -> 1 
Repeat it, to file: setup. h 

4 At system variable [path] : 

[My Computer] [Properties] [Advanced] 

[Environment Variables] [System Variables] 

Path [Edit ] [Variable value], Add: 

" ; X : \Program Files\CodeBlocks\MinGW\bin; " 

5 Command Prompt : >cd\ 

>cd X : \wxWidgets-2 . 8 . ll\build\msw 
>mingw32-make -f makefile. gcc MONOLITHIC=0 
S HARED =0 UNICODES USE_OPENGL=l 
BUILD=release 

Re-Command, with option: 

BUILD=release -> BUILD=debug 

6 Download (~100KB) , Extract (~700KB) : 
http : //code . google . com/p/rga/ downloads/ 

VEMSA3D_source_l . zip -> X : \ 

7 IDE CodeBlocks, [File] [Open] : 

X : \VEMSA3D\build\win-cb\VEMSA3D . workspace 

* If system disk is not C:\ then: 

(#) [Project] [Build options. . .] [VEMSA3D] 

* [Search directories] Correct in each of: 

* [Compiler] , [Linker] , [Resource compiler] 

* C:\ -> X : \ [Yes] 

* Repeat from (#) , for [Release] 

* Repeat from (#) , for [Debug] 

[Build] [Select target] [Release] [Build] 

[Build] [Select target] [Debug] [Build] 
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This process results in 
[VEMSA3D] executables of ~35 
MB Debug and of ~5MB Release 
versions, which run under NT4, 
W2K and WXP, at least. 

Alternatively, the use of 
the freeware IDE MS Visual 
Studio Express with C++ 2008 
Compiler is also suggested. 
For that, first download MS 
Visual Studio 2008 Express 
iso-image from Microsoft web- 
site, burn it into a CD-R, 
and install it. Then, down- 
load from our repositories 
the file [VEMSA3Dfiles-4win- 
cb.zip], extract it, and fol- 
low the included setup in- 
structions. This process re- 
sults in [VEMSA3D] executa- 
bles of -5MB Debug and of 
~2MB Release versions, which 
run under W2K and WXP, at 
least, but definitely do not 
run under NT4 . 

Practical Antenna Applica- 
tions and Results 

This section presents the 
implementation of [VEMSA3D] to 
produce eight antenna models, 
from simple to more compli- 
cated structures. The input 
data were imported through 
[Import Geometry] of the [Fi- 
le] menu (Fig. 2.1). 

Fig. 7 shows an array of 2 
dipoles for operation at the 
frequency of 1111 MHz, dis- 
tanced by 0.85A, constructed 
by bare copper wire of 1 mm 
(0.0037A) radius, and mea- 
sured. In the same figure, 



the 3D radiation intensity 
pattern and its 2D main plane 
cuts are shown [11] . The con- 
tinuous line represents [VEM- 
SA3D] results, the dashed li- 
ne, analytically produced pat- 
terns, and the dots, measure- 
ments made using the authors' 
group VNA system [12], [13]. 

Fig. 8 illustrates the re- 
sults for a constructed and 
measured improved Hentenna mo- 
del, at 1110 MHz, with height 
of A/2, width of A/6, and 
with the active element at a 
distance 7A/12 from its bot- 
tom [14] . 

The third antenna consists 
of a A/4 monopole (5.83 cm) 
at 1286 MHz over a circular 
counterpoise of 14 cm radius. 
The monopole was constructed 
by 2 mm diameter copper wire 
and the counterpoise has been 
printed on a circuit board of 
29 cm x 29 cm, as a circle 
with four radials of 2 mm 
width [15]. In Fig. 9, the 
left 3D and 2D patterns are 
for vertical polarization, 
while the right ones for ho- 
rizontal. For a better repre- 
sentation, the patterns have 
been normalized. 

A somehow more complicated 
antenna is shown in Fig. 10. 
The left 3D and 2D patterns 
correspond to the vertical po- 
larization of a monopole over 
circular counterpoise with 16 
radials, while the right co- 
lumn patterns are the corre- 
sponding 3D and 2D patterns 
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for the same antenna, but with 
an additional grounded disc 
below it, and without elec- 
trical connection to the coun- 
terpoise antenna, at a dis- 
tance of 0.9 cm. The grounded 



disc was constructed from thin 
copper sheet and was modeled 
in [VEMSA3D] with 5 concen- 
tric circles and 64 radials 
[15]. 




Fig. 7: Array of 2 dipoles: 73 dipole current modes 




— VEMSA3D •••Measurements 



Fig. 8: Hentenna: 109 dipole current modes 
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Fig. 9: Monopole over a counterpoise with 4 radials: 
167 dipole current modes 




— VEMSA3D •••Measurements 



Fig. 10: Monopole over a counterpoise with 16 radials: 
419 dipole current modes. 

Monopole over a counterpoise with 16 radials over 
a grounded disc: 2594 dipole current modes 
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In Fig. 11, the results of 
[VEMSA3D] for a commercial 
UHF antenna are shown [16] . 
The antenna model is presen- 
ted in detail and separately, 
for the active element sys- 
tem, as well for the wire- 
frame reflector. The patterns 
correspond to the center fre- 
quency of its operation at 
650 MHz. 

Fig. 12 illustrates the 
model of a small jet airplane 
from the well-known freeware 
[4NEC2X] antenna simulator 
[17] . Simulation was carried 
out at the frequency of 10 
MHz with the same number of 
points and segments in both 
simulators. There is a good 
agreement between them for 



the radiation intensity pat- 
terns although some deviation 
in input impedance and direc- 
tivity is observed. 

Fig. 13 shows the results 
for a horn antenna with di- 
mensions proposed by K. Pitra 
and Z. Raida [18]. A small 
bow-tie feeder with a trian- 
gular perimeter of 0.47A and 
flare angle 37.50°, is con- 
sidered. The antenna was ini- 
tially simulated at the fre- 
quency of 40 GHz in both 
[4NEC2X] and [VEMSA3D] simu- 
lators, with 1205 segments to 
be consistent with the re- 
stricted maximum number of 
1500 segments of [4NEC2X] . 
There is a good agreement be- 
tween the produced results. 




Fig. 11: Typical commercial 
723 dipole current 



TV UHF antenna: 
modes . 
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[VEMSA3D] 



Fig. 12: Jet airplane modeled in [4NEC2X] and [VEMSA3D] : 
391 dipole current modes 



Finally, the most compli- 
cated antenna model, which is 
presented, corresponds to the 
same horn antenna at the same 
frequency of 40 GHz, with a 
dense wire-frame consisting 
of 3266 points and 4362 seg- 
ments. The model is shown in 
Fig. 14. The total number of 
dipole current modes is 5458, 



as shown in Fig. 15. The pro- 
cess time ranges from ~40 min 
in an AMD Phenom X2 550 3.11 
GHz CPU to -220 min, in an 
Intel Pentium 4 1.7 GHz CPU. 
Fig. 16 illustrates the re- 
sulting radiation patterns 
and Fig 17 the current ampli- 
tude and phase. 
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Ji 

w 

[4NEC2X] 




[VEMSA3D] 




Impedance | 


30 -i 


0.06 


S.W.R. 50 j 


1.67 




Frequency 


1 40000 


Mhz 


Wavelength 


| 7.e-3 


mtr 


Norm-All : 


: 6.53 dBi 






Resistance [Ohm] 

4 . 673E+001 

Reactance [Ohm] 

-3 . 929E+000 

VSWR [] 

1.1115 

Directivity 
6. 702E+000 [dBi] 



Fig. 13: Horn antenna with bow-tie feeder modeled in 
[4NEC2X] and [VEMSA3D] : 1370 dipole current modes 




Fig. 14: Horn antenna with bow-tie feeder: 4362 segments 
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Points Segments Dipole Current Modes ; input Sources Power 





Point 1 


Point 2 


Point 3 


Segment A 


Segment B 


Amplitude [A] 


Phase [deg.] 


5455 


1121 


3263 


1122 


4319 


4320 


1 . 053E-002 


-170.85 


5456 


1121 


3264 


1123 


4321 


4322 


1 . 042E-002 


-172.78 




1122 


3265 


1123 


4323 


4324 


1 . 615E-004 


-127.73 


^5458 


) 1121 


3266 


1118 


4325 


4326 


2 . 602E-002 


12.29 



Points Segments Dipole Current Modes 



Total Input Power (Real) 


2 . 543E-002 [W] 


Total Input Power (Imag) 


-5 . 540E-003 [W] 


Total Radiated Power 


2 . 543E-002 [W] 


Total Dissipated Power 


0 . 000E+000 [W] 


Power Efficiency 


100% 


Power Gain 


5 . 719E+000 [] 




7 . 573E+000 [dBi] 




5 . 424E+000 [dBd] 


Directivity 


5.719E+000 [] 




7 . 573E+000 [dBi] 




5 . 424E+000 [dBd] 



Fig. 15: Horn antenna with bow-tie feeder: 
5458 dipole current modes 




Fig. 16: Horn antenna with bow-tie feeder: 
Radiation intensity 
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Fig. 17: Horn antenna with bow-tie feeder: 
Current Amplitude and Phase 



Conclusion 

The first stable version 
of a visual EM simulator has 
been developed and released 
as FLOSS. The computational 
results of its use were found 
to be in good agreement with 
experimental measurements as 
well with the comparable free- 



ware simulator [4NEC2X] . Sin- 
ce the number of possible mo- 
difications and additions to 
the attributes of [VEMSA3D] 
seems to be endless, no at- 
tempt will be made to suggest 
a particular direction for its 
future development: any con- 
tribution from the antenna 
community is very welcomed. 
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Abstract 

The aim of this paper is the study of various parameters 
affecting the lifetime of transformer oil in transformers of 
150/20 kV. Fifty (50) samples of oil were taken from such 
transformers in the major Athens area, Greece. The parame- 
ters investigated - according to international standards - 
were breakdown strength, oil color, humidity, interfacial 
tension and tan5. Thermal and mechanical stresses have as 
result the oxidation of transformer oil and the deteriora- 
tion of its insulating properties. Humidity and foreign par- 
ticles also consist factors contributing to the lowering of 
the breakdown strength of transformer oil. In most examined 
samples the breakdown strength and tan5 were satisfactory. 
It is shown that the slightest contact with the atmospheric 
air may affect humidity. A color index of 3 does not neces- 
sarily mean that the oil is bad. Foreign particle presence 
combined with humidity may decrease the interfacial tension. 
Generally speaking, the 50 investigated transformers were in 
a satisfactory state and none of them was required to have 
its oil replaced. A main conclusion of this work is that we 
should not base our judgment about the oil quality on only 
one or two parameters but on a combination of more parame- 
ters . 
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Introduction 

Transformer oil is a very 
important component of a 
transformer. It must have 
good thermal and insulating 
properties [1] . The insulat- 
ing oil is used for providing 
insulation between the live 
parts of the transformer and 
the grounded parts as well as 
for carrying out the heat 
from the transformer to the 
atmosphere [2] . Faults in 
transformers are rare (l%-2% 
per year in power transmis- 
sion and distribution sys- 
tems), but they can be very 
costly in terms of economic 
and technical consequences. 
Faults can due, among others, 
to lightning and switching 
overvoltages, insulation fai- 
lure, humidity, foreign par- 
ticles and bad maintenance. 
The latter is a significant 
factor for the lifetime of a 
transformer [3], [4]. The main 
factors which may accelerate 
the ageing of the oil are hu- 
midity, temperature varia- 
tions, oxidation and foreign 
particles. Various factors af- 
fecting the breakdown strength 
of transformer oils have been 
discussed and analyzed [5], 
[6] . There is no single mea- 
surement which can deliver 
enough information as to the 



ageing and/or deterioration of 
transformer oil, mainly be- 
cause transformer insulation 
is a dynamic system, in which 
e.g. humidity may migrate 
from the oil to paper insula- 
tion and from paper insula- 
tion back to the oil [3], 
[7] . With this in mind, a va- 
riety of diagnostic methods 
were employed in order to stu- 
dy the state of fifty trans- 
formers of 150/20 kV. The who- 
le work was carried out with 
the aid of the Public Power 
Corporation (PPC) Transformer 
Division in Athens, Greece. 
The fifty investigated trans- 
formers were from the major 
Athens area. 

Diagnostic Methods 

Warning signs about the 
state of a transformer are, 
among others, a big increase 
of partial discharges (» 
2500 pC), a visible deterio- 
ration because of foreign me- 
tallic and carbon particles, 
the presence of humidity in 
the solid insulation about 3- 
4% and the presence of slud- 
ge. The latter is the last 
visible state of deteriora- 
tion. Experience indicates 
that the breakdown behavior 
and breakdown voltage are de- 
termined much more from the 
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above mentioned factors than 
that of pure insulating li- 
quids [4], [8], [9] . 

Several diagnostic methods 
were used in order to see the 
quality of the transformers 
in question. The characteri- 
zation of the oil color (DIN 
51517 - ASTM 155) was per- 
formed through a device (chro- 
mometer) including standard 
glass disks and two glass 
jars with lid. The control of 



breakdown voltage was mea- 
sured by a typical Foster 
test cell, according to IEC 
156/95 (Fig. 1). 

The control of humidity in 
the oil was measured by a 
Metrohm - 684 KF Coulometer, 
which consisted of a glass 
container with a stirrer ti- 
tration in which the reagent 
from container storage is ad- 
ded . 




Fig. 1: The test cell for breakdown voltage measurements 
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The device is fully auto- 
mated and once the experi- 
menter gives the settings, it 
measures the moisture content 
of the oil. The measurements 
were performed according to 
IEC 814. The control of in- 
terfacial tension (ASTM D971 
- 91) was performed via a 

tensimeter, which gives the 
value in dynes per centimeter 
in a direct reading. 

The device that performed 
measurements of tan5 and of 



resistivity, is the BAUR-DTL 
fully automated device for 
measuring dielectric losses 
of oils. Such a system has a 
fully automated process for 
measuring dielectric loss, re- 
lative dielectric constant 
and resistivity (Fig. 2). The 
measurements were performed 
according to IEC 247. The 
density of oil was performed 
according to DIN 51517, with 
the aid of a pipe of 250 ml, 
an electronic thermometer and 
a glass cylinder. 




Fig. 2: Device for measurement of tan5 and resistivity 
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It is true that no single 
diagnostic method can give 
full information as to the 
state of the transformer oil. 
The aforementioned methods 
may give a better picture of 
its state. 

Results 

The sampling was done ac- 
cording to specification ASTM 
D 923. Sampling should take 
place in clean conditions (ab- 
sence of humidity and pollu- 
tion), suitable glass vessels 



should be used and the latter 
should be kept clean and her- 
metically closed. Every sam- 
ple should be kept away from 
light according to VDE 0370/ 
9.61. Every sample should be 
taken while the oil is warm. 

In Tab. 1 a classification 
of values of the various in- 
vestigated parameters of in- 
sulating oil is given. An oil 
can be classified as good, 
acceptable or bad according 
to Tab . 1 [10] . 



Tab. 1: Classification of insulating oils 



Oil Parameters 


Good 


Acceptable 


Poor 


Color 


< 2 


- 


> 2 


Breakdown Voltage (kV) 


> 40 


30 - 40 


< 30 


Humidity [ppm] 


< 10 


10 - 25 


> 25 


Interfacial Tension (mN/m) 


> 28 


22 - 28 


< 22 


tan5 


< 0.1 


0.1 - 0.5 


> 0.5 


Resistivity (p) (GQ-m) 


> 3 


0.2 - 3.0 


< 0.2 



The sampling of oil from 
the fifty transformers 150/20 
kV was done with due care and 
according to the standard 
practice. In Figures 3-7, the 
results of the measurements 
are shown, regarding tan5, co- 
lor, humidity and interfacial 



tension respectively. Green 
color symbolizes the good sam- 
ples, yellow color the ac- 
ceptable samples whereas the 
red color shows the bad sam- 
ples. From Figs. 3 and 4, it 
is clear that the breakdown 
voltage and tan5 values of 
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most of the samples are very 
good. This is due to syste- 
matic control of the oil and 
the good maintenance. Figs. 3 
and 4 indicate that these two 
factors, which are related to 
ageing and oxidation, are 
relatively stable. 

In Fig. 5, humidity is in 
relatively acceptable levels. 
Only a small percentage of 
the transformers (4%) seems to 
have high humidity. Most of 
the samples are within the 
limits prescribed by the in- 
ternational standards. It 
must be emphasized that the 
humidity level is a parameter 
which changes easily, since 
the slightest contact of the 
oil with the atmosphere may 
change its characteristics. 
In Figs. 6 and 7 the results 
regarding the color as well 
as the interfacial tension 
are shown. 

As time passes by, the oxi- 
dation products change the 
oil color. Most of the sam- 
ples had a rather acceptable 
color. Even a color index of 
3 does not necessarily con- 
sist an objective indication 
of the oil quality. For this 
reason, color measurements 
should in fact be accompanied 
by other parameter measure- 
ments. In fact, although 36% 
of the samples showed a ra- 
ther dark color (Fig. 6) other 
parameter measurements indi- 
cated that these samples were 
good or acceptable. The pre- 



sence of foreign particles in 
combination with humidity may 
reduce the interfacial ten- 
sion of the oil. In the in- 
vestigated samples, a per- 
centage of 28% (Fig. 7) is 
characterized as poor. This, 
however, is not particularly 
annoying, if we take into ac- 
count for these samples also 
the other parameter measure- 
ments . 

The density of the inves- 
tigated oil samples was mea- 
sured in the generally ac- 
ceptable values, i.e. between 
0.85 and 0.92 gr/ml (with the 
lowest recorded being 0.85 
gr/ml, whereas the highest was 
0.91 gr/ml). Although the oil 
density does not consist per 
se an individual characteris- 
tic of the examined sample, 
its increase may imply an in- 
crease of degradation by- 
products. In the context of 
the present work the oil den- 
sity was used for the calcu- 
lation of the interfacial 
tension [11] , [12] . 

It can be said that, in 
general, the state of the in- 
vestigated transformer oil 
samples was more or less sa- 
tisfactory. In a few cases, 
there is a need of further 
filtering and possibly a se- 
cond sampling it should be 
carried out. Whereas no trans- 
former functioned with a par- 
ticularly bad oil, it is true 
that, with transformer age- 
ing, the oil suffers from so- 
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lid impurities, free and dis- 
solved water particles and 
dissolved air. Frequent sam- 
pling is necessary in order 
to ensure the good function- 
ing of such transformers. The 
results reported here are in 
line with those published be- 
fore [11] . Although a statis- 
tical approach of the whole 
subject is desirable (i.e. to 
try to correlate the data 
collected here as well as 
from previous published work 
with the pre-history of each 
individual transformer), it 
is difficult to be realized 
since the Greek Electricity 
System has transformers from 
a variety of manufacturers. 



This inhomogeneity of suppli- 
ers certainly renders the re- 
lation between the data col- 
lected with the pre-history 
of transformers very diffi- 
cult. It also should be noted 
that with the term "pre-his- 
tory" we mean the detailed 
registration of all faults, 
faulty conditions, lightning 
strokes, switching overvolt- 
ages etc. which have occurred 
in a transformer. Previous 
work done in this direction 
was only partially successful 
since at that time, pre-his- 
tory of transformers was only 
related to one parameter, 
namely that of breakdown 
strength [13], [14], [15]. 



2% 




Fig. 3: Graph of breakdown voltage results 
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2% 





Fig. 5: Graph of humidity results 
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It can be said that this 
piece of work here does not 
consist per se an original 
piece of research. This paper 
does not claim to have ex- 
plored new inroads regarding 
the mechanisms of breakdown 
of dielectric liquids, as for 
example in [16], or novel di- 
rections regarding new diag- 
nostic methods, as for exam- 
ple in [17]. The whole pur- 
pose of this paper is, how- 
ever, to show that monitoring 
work is necessary - not to 
say essential - for the cor- 
rect maintenance of trans- 
former oil in bigger trans- 
formers. An attempt to also 
relate the present state of 
transformer oil in the above 
mentioned transformers with 
their previous history will 
follow. 

Some Further Remarks 

One may ask which from the 
parameters investigated, are 
the most important. It is a 
difficult question to answer. 
However, if one has to choose 
between the aforementioned pa- 
rameters, he would most pro- 
bably select two of them, na- 
mely humidity and breakdown 
strength. If one looks care- 
fully Fig. 3 (referring to 
breakdown voltage results) 
and Fig. 5 (regarding humid- 
ity results), one may see 
that 98% of the investigated 
transformers had good oil 
breakdown strength and 96% of 



the investigated transformers 
had good or acceptable levels 
of humidity. This means that 
Fig. 3 and Fig. 5 are in more 
or less good agreement. Hu- 
midity plays a critical role, 
since it can contribute to a 
dramatic lowering of break- 
down strength, as was also in- 
dicated in some older but 
nevertheless classical publi- 
cations [18], [19], [20]. Ge- 
nerally speaking, inclusions 
of humidity more than about 

10 ppm (at normal tempera- 
ture) cause a lowering of 
breakdown strength. A low oil 
breakdown strength may imply 
that there are foreign parti- 
cles and/or admixtures in the 
oil. On the other hand, a 
high breakdown strength does 
not necessarily mean that the 

011 is good. It may be possi- 
ble that the quantity of fo- 
reign particles may not be 
sufficiently large, so that 
it can influence the break- 
down voltage [21] . 

The interfacial tension, 
although it gives an idea of 
the concentration of oxida- 
tion by-products in a trans- 
former oil, is not necessa- 
rily an indicator for defi- 
nite conclusions about the 
oil under investigation. This 
is because in warmer periods 
of the year, the oil tempera- 
ture increases and humidity 
may affect the oil more than 
in cooler periods [22] . 

Change of the oil color 
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may mean the existence of by- 
products or the presence of 
foreign particles. Although 
such a change may imply a 
certain degree of pollution 
of the oil under investiga- 
tion, the color by itself may 
not be considered as a very 
reliable indicator of the oil 
quality, as was shown in this 
paper and was also reported 
in [23]. 

Tan5 results (Fig. 4) 
match very well with Fig. 3 
results of breakdown 
strength. Tan5 changes as the 
oil degrades. Although tan5 
measurements cannot be taken 
as a sole criterion of oil 
quality, its results match ex- 
tremely well with the break- 
down strength data, in the 
context of the present paper. 
The low value of tan5 depends 
on the nature of the oil as 
well as on its processing 
[22], [24], 

A last remark should be 
made concerning the monito- 
ring of both transformers of 
150/20 kV and of distribution 
transformers (both kinds of 



transformers in the major 
Athens area) : it seems that 
the former have a larger per- 
centage of good/acceptable oil 
than the latter [25] . To va- 
lidate, however, this point, 
more work has to be done with 
transformers of both kinds 
from the major Athens area. 

Conclusion 

Insulating oil samples 
from transformers of 150/20 
kV have been investigated. In 
the context of this work, se- 
veral parameters - through 
the appropriate diagnostic 
techniques - affecting the 
state and lifetime of trans- 
former oil have been studied. 
No single parameter can fully 
describe the state of the oil 
of a transformer. The variety 
of parameters investigated he- 
re may give a more complete 
picture. In the context of 
this work, the majority of 
the investigated transformer 
oil samples were found to be 
good or acceptable. This 
points out to the further con- 
tinuing sampling at regular 
intervals . 
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